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In this paper we derive by direct calculation the statistical properties of extremes of general 
observables for chaotic system, taking as mathematical framework that given by Axiom A dynamical 
systems possessing an invariant SRB measure. We prove that the extremes of so-called physical 
observables are distributed according to the classical Generalised Pareto Distribution (GPD) and 
derive explicit expressions for the scaling parameter a and the shape parameter £. In particular, 
we derive that £ has a universal expression which does not depend on the chosen observables, with 
£ = — l/(d s + d u /2 + d„/2), where d s , d u , and d„ are the partial dimensions of the attractor of the 
set in the stable, unstable, and neutral directions, respectively. The shape parameter is negative and 
is close to zero when high-dimensional systems are considered. This result agrees with what derived 
recently using the Generalised Extreme Value approach. Combining the results obtained using such 
physical observable and the properties of the extremes of distance observables considered in previous 
papers, where a direct link is found between the £ parameter and the Kaplan- Yorke dimension of 
the attractor, it is possible to derive estimates of the partial dimensions of the attractor along 
the stable and the unstable directions of the flow. Moreover, by rewriting the expression of the 
shape parameter as suitable combination of moments of the considered observable and by using 
the Ruelle response theory, we provide a general framework for describing the sensitivity of the 
shape parameter with respect to e-perturbations to the flow. The possibility of defining the linear 
response for £ with respect to e provides mathematical arguments for suggesting that the Kaplan- 
Yorke dimension of the attractor changes with a certain degree of smoothness when small external 
perturbations are introduced, as typically found in numerical experiments with intermediate to 
high-dimensional chaotic systems. 



I. INTRODUCTION 

Extreme value theory (EVT) is gaining more and more 
prominence in a vast range of scientific fields because 
of its theoretical relevance in mathematical and physical 
sciences, and because it addresses the problem of under- 
standing, modeling, evaluating risk factors such as those 
related to instabilities in the financial markets and to nat- 
ural hazards related to seismic, climatic and hydrological 
extreme events. Even if the probability of extreme events 
is very low and decreases quickly with their magnitude, 
the associated risks can dominate over those coming from 
events belonging to the bulk of the statistics. An exten- 
sive account of recent results and relevant applications is 
given in [1]. 

The first comprehensive discussion of EVT dates back 
to the fundamental work by Gnedenko [2], who inves- 
tigated the distribution of the maxima of a sample of 
independent identically distributed (i.i.d) stochastic vari- 
ables. He showed that under very general conditions such 
maxima are distributed according to the so-called Gen- 
eralised Extreme Value (GEV) distribution. The clas- 
sic way of dealing with the statistical inference of ex- 
tremes actually follows quite precisely the steps of the 
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Gnedenko's theorem. One partitions the experimental 
time series into bins of fixed length, extracts the max- 
imum of each bin, and fit these data to the GEV dis- 
tribution family using, typically, methods such as max- 
imum likelihood estimation (MLE) or L-moments. See 
[3] for a detailed account of this methodology. It is pos- 
sible to deal with extremes by taking a different point of 
view, i.e., by defining extremes as the events exceeding 
a given threshold. In the limit of very high threshold, 
we expect that the extremes are distributed according to 
the Generalized Pareto Distribution (GPD) introduced 
by Pickands III [4] and Balkcma and De Haan [5]. In 
the case of i.i.d. variables, it is well known that a strong 
connections exists between the two methodologies. As 
shown in [6] , we have that if block maxima obey the GEV 
distribution, then exceedances over some high threshold 
obey an associated GPD. Nonetheless, it is apparent that, 
whereas the two approaches are equivalent in the asymp- 
totic limit, the GPD approach provides more reliable and 
more robust results when realistic, finite time scries arc 
considered (see, e.g., [7]). 



A. A Brief Recapitulation of Extreme Value 
Theory 

First, we briefly recapitulate the GEV approach to 
EVT. Gnedenko [2] studied the convergence of maxima 
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of i.i.d. variables 



Xn, Xi , ..X, 



m—l 



with cumulative distribution function (cdf) F m (x) — 
P{a m (M m — b m ) < x} where a m and b m are normaliz- 
ing sequences and M m = max{A , Xi, X m _i}. Under 
general hypothesis on the nature of the parent distribu- 
tion of data, Gnedenko [2] showed that the asymptotic 
distribution of maxima, up to an affine change of vari- 
able, belongs to a single family of generalized distribu- 
tion called GEV distribution whose cdf can be written 



as: 



where 



lim F m (x) = F GEV (x;fi,a,K) = e t{x) (1) 



This expression holds for 1 + k(x — fi)/a > 0, using (igl 
(location parameter) and a > (scale parameter) as scal- 
ing constants, and k e R is the shape parameter (also 
called the tail index). When k — )■ 0, the distribution cor- 
responds to a Gumbel type (Type 1 distribution) . When 
the index is positive, it corresponds to a Frechet (Type 2 
distribution); when the index is negative, it corresponds 
to a Weibull (Type 3 distribution). 

We briefly mention the Pareto approach to EVT. We 
define an exceedance as z = X — T, which measures by 
how much the variable X exceeds a given threshold T. As 
discussed above, under the same conditions under which 
the block maxima of the i.i.d. stochastic variables X 
obey the GEV statistics, the exceedances z are asymp- 
totically distributed according to the Generalised Pareto 
Distribution. Defining F T (z) = P(X - T < z\X > T), 
we have that liniT^oo Ft(z) — Fqpd{z; £, cr), with [4, 6]: 



F GPD (z;£,a) 



l-exp(-f) 



-i/e 



for £ ? 0, 
for £ = 0, 



(3) 



where the range ofzis0<z<ooif£>0 and < z < 

if t < o. 

The relation between GEV and GPD parameters has 
been already discussed in literature in case of i.i.d vari- 
ables [7-10]. It is first interesting to note that 

F GPD (z-T;a,0 = 1 + log (F GEV (z; T, a, £)) (4) 

for F GEV (z; fi, a, k) > exp -1 , where the latter condition 
implies z > T [11]. If we consider the upper range z 3> T, 
we have that F GE v(z; T,a,£) is only slightly smaller than 
1, so that Eq. 4 implies that 

F GPD (z - T;a,C) - F GEV (z;T,a,£,) 

+ 0((l + £(z-T)A7)- 2/ «), 



so that the two distributions are asymptotically equiv- 
alent. This simple result is actually equivalent to the 
rather cumbersome formulas given in Coles [9] and Katz 
et al. [8] (and reported also by us in [12] for defining the 
correspondence between the parameters of the GPD and 
GEV distributions describing the statistics of extreme 
events extracted from the same data series. 



B. Extreme Value Theory for Dynamical Systems 

Recently, a great deal of attention has focused on un- 
derstanding to what extent EVT can be applied to study 
the extreme of observablcs of deterministic dynamical 
systems. The main applications-driven motivation for 
this renewed interest probably comes from the spectacu- 
lar development of numerical modeling in a geophysical 
fluid dynamical context and from the need to assess the 
ability of climate model to reproduce the observed statis- 
tical properties of extremes in present climate conditions 
and understand how they will change in an altered cli- 
mate [13]. Other related applications include the numer- 
ical simulation of hydrological risk and of the production 
of electric energy from wind. It is clear that the matter is 
far from being trivial: numerical experiments on climate 
models of various degrees of complexity have shown that 
the speed of convergence (if any) of the statistical prop- 
erties of the extremes definitely depends strongly on the 
chosen climatic variable of interest [3, 14-16]. 

Apart from these specific albeit very relevant applica- 
tions, this problem has been addressed by the mathemat- 
ical and statistical physical community. A first important 
result is that when a dynamical system has a regular (pe- 
riodic of quasi-periodic) behaviour, we do not expect, in 
general, to find convergence to GEV distributions for the 
extremes of any observable Nicolis et al. [17], Haiman 
[18]. Instead, if one choses specific observables and con- 
siders dynamical systems obeying suitable mixing condi- 
tions, which guarantee the independence of subsequent 
maxima, it is possible to prove that the distribution of 
the block maxima of the observables converge to a mem- 
ber of the GEV family. The observables are expressed as 
g(dist(x,Xo)), a function g of the distance of the orbit x 
from a point in the attractor xq, usually taken as the ini- 
tial condition, such that g(y) has a global maximum for 
y = 0. The specific member of the GEV family (which 
is determined by the sign of the shape parameter) the 
maxima distribution converges to depends on the specific 
choice of g. The paper by Collet [19] can be considered 
the cornerstone for the subsequent results obtained in the 
last few years [20-22] . The resulting parameters of the 
GEV distributions can be expressed as simple functions 
of the local dimension of the attractor. These results have 
been shown to be accurately detectable in numerical ex- 
periments when considering finite time series [23-25] . If, 
instead, the maxima are clustered, so that they feature 
a relative strong short-time correlation, the results have 
to be modified by introducing the extremal index [26]. 
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Recently, it has been shown how to obtain results 
which arc independent on whether the underlying dy- 
namics of the system is mixing or, instead, regular . The 
key ingredient relies on using the Pareto rather than the 
Gnedenko approach. Such a shift in the point of view on 
extremes allows to derive results that do not dependent 
on whether extremes feature strong time-correlations or 
not. Assuming only that the local measures scales with 
the local dimension [27], it is possible to obtain by direct 
integration a GPD for the threshold exceedences of the 
observables g(dist(x,x )) introduced in [20-22, 26] when 
considering a generic orbit of a dynamical systems. In 
fact, the Pareto approach entails sampling all the points 
of the orbit that are very close to xq, thus sampling the 
local scaling of the invariant measure. With a suitable 
choice of g(dist(x,Xo)), the resulting £ of the GPD is 
proportional to the inverse of local dimension [12]. The 
results obtained using the Pareto approach agree exactly 
with what derived using the Gnedenko approach under 
the assumption of mixing dynamics in [20-25]. When 
the underlying system is mixing enough, the dynamical 
(Gnedenko approach) and the geometrical (Pareto ap- 
proach) points of view on extremes give the same results, 
whereas differences emerge if the dynamics is such that 
strong time correlations exist between block-maxima. 
The selection of the extremes of the g observables dis- 
cussed above - using either the Gnedenko or the Pareto 
approaches - acts as magnifying lens near the initial con- 
dition, and that's why one can extract information on the 
local dimension. Therefore, this provides a potentially 
viable alternative to e.g. the Grassberger-Procaccia al- 
gorithm [28] for the investigation of the scaling properties 
of the invariant measure of a chaotic attractor. 

The results discussed above feature a major drawback 
when considering their relevance in many applications. 
Extreme events corresponds to close returns of the orbit 
to to its initial condition. While relevant problems in 
natural sciences can be set in the framework of this class 
of observables (e.g. the classic problem of 'weather ana- 
logues" in meteorology, already discussed by Lorenz in 
connection to the problem of predictability [29]), this is 
not the typical structure of the observables encountered 
in many applications, such as the case of total energy 
or enstrophy of a fluid flow. Recently, this problem has 
been addressed in Holland et al. [30], who have stud- 
ied, using the GEV approach, whether EVT applies for 
physical observables of maps obeying the mixing condi- 
tions proposed in Freitas and Freitas [20], Freitas et al. 
[21], Gupta et al. [22], Freitas et al. [26]. They consider a 
general observable A = A(x) reaching its maximum value 
A mo .x restricted to the support of the invariant measure 
in x — x (assuming for simplicity that such a point is 
unique), and assume that VAI^^ ^ 0. We note that 
A(x) is indeed not of the form g(dist(x, xo)) discussed 
above. They find that the block maxima of A(x) are 
asymptotically distributed according to a member of the 
GEV distributions, where the shape parameter is nega- 
tive and can be written as a simple function of the par- 



tial dimensions along the stable and unstable manifolds 
at the point. This seems indeed a very relevant result, as 
it provides a universal property of extremes, regardless 
of the specific functional form of the observable A. 



C. Goals of the Paper 

In this paper, we consider a GPD approach to EVT 
and try to complement and improve the results presented 
in [30] regarding the physical observables and those pre- 
sented in [12] regarding the distance observables. We 
focus our attention on Axiom A systems [31], which are 
a special class of dynamical systems possessing an SRB 
invariant measure [32] and featuring hyperbolicity in the 
attracting set. Such invariant measure coincides with the 
Kolmogorov's physical measure, i.e. it is robust against 
infinitesimal stochastic perturbations. Another impor- 
tant property of Axiom A systems is that it is possible to 
develop a response theory for computing the change in 
the statistical properties of any observable due to small 
perturbations to the flow [33, 34]. Such response the- 
ory has recently been the subject of intense theoreti- 
cal [35, 36], algorithmic [37] and numerical investigations 
[38-41] and is gaining prominence especially for geophys- 
ical fluid dynamical applications. Moreover, the response 
theory seems to provide powerful tools for studying mul- 
tiscale systems and deriving parametrizations of the im- 
pact of the fast variables on dynamics of the slow [42, 43]. 
Finally, an important property of Axiom A systems is 
that, while the dynamics of natural or artificial systems is 
definitely not Axiom A in general, Axiom A systems can 
be considered as good 'effective' models of actual systems 
with many degrees of freedom thanks to the so-called 
chaotic hypothesis, which is somewhat the equivalent in 
the non-equilibrium framework of the classic ergodic hy- 
pothesis for equilibrium dynamics [44] . Moreover, as dis- 
cussed in [41], when we perform numerical simulations we 
implicitly assume that the system under investigation is 
Axiom A or Axiom A-equivalent. Therefore, considering 
Axiom A systems seems a good mathematical framework 
in view of providing results useful for a large spectrum 
of applications. The choice of considering Axiom A sys- 
tems is instrumental in the derivation of various results 
on the relationship of EVT parameters to the dynamical 
and geometrical properties of the system, and will allow 
addressing the problem of the sensitivity of extremes to 
small perturbations of the system. The dependence of the 
properties of extremes of parametric modulations of the 
underlying dynamics is an issue of relevant theoretical as 
well as applicative interest. The practical interest stems 
from the fact that it is relevant to be able to control or 
predict variations in extreme events due to small pertur- 
bations to the dynamics. The theoretical interest comes 
from the fact that when considering extremes, univer- 
sal parametric probability distributions can be defined, 
as opposed to the case of the bulk statistical properties. 
Because of this, we may hope to reconstruct the parame- 
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ters descriptive of the EVT from simple moments of the 
distributions, express these in terms of observables of the 
system, and use the Ruelle response theory for express 
rigorously the sensitivity of extremes to small perturba- 
tions to the dynamics. 

In Sec. II we show that by direct integration it is pos- 
sible to derive the value of the two GPD parameters £ 
and a, and, in particular, that the value of £ agrees with 
the GEV shape parameter obtained in [30] . We also show 
that, combining the results obtained using such physical 
observable A and the distance observables considered in 
[12], it is possible to derive the estimates of the partial 
dimensions of the attractor along the stable and the un- 
stable directions of the flow. In Sec. Ill, we develop a 
linear response theory describing the impact of small e- 
perturbations to the flow on the statistical properties of 
the extremes of the observable A. In particular, we will 
focus on computing the changes of the shape parameter £. 
We will find two equivalent expression for the sensitivity 
of £ with respect to e. First, we will provide an expression 
for the sensitivity of £ in terms of the first two moments 
of the probability distribution of above-threshold A(x) 
occurrences. Such expression entails a combination of ob- 
servables of the Axiom A system, so that one can use Ru- 
elle's theory to compute the response to e-perturbations 
to the dynamics. Nonetheless, we show that mathemat- 
ical problems emerge when considering some limits. We 
follow the same approach for defining an expression for 
the sensitivity of the extremes of the distance observables 
considered in [12]. Then, we will relate the sensitivity of 
£ to the sensitivity of the Kaplan- Yorke dimension of the 
attractor. We will link our results to the del-established 
fact that both the extremes of observables and quantities 
like the Lyapunov exponents feature a good degree of reg- 
ularity with respect to perturbations when one considers 
intermediate complexity to high-dimensional chaotic dy- 
namical systems. In Sec. IV we comment the obtained 
results and present our conclusions. 



II. EXTREME VALUE THEORY FOR 
PHYSICAL OBSERVABLES OF AXIOM A 
SYSTEMS 

A. Geometry of the problem 

Let us consider a continuous-time mixing Axiom A 
dynamical system x = G{x) on a compact manifold 
M C M d , where x{t) = f{x in ), with x(t = 0) = x m e M 
initial condition and /* evolution operator, is defined for 
all t € M>o- Let us define f2 as the attracting invari- 
ant set of the dynamical system , so that v is the asso- 
ciated SRB measure with support O = supp(V). Let 
us now consider a smooth observable A whose maxi- 
mum restricted to the support of v is unique, so that 
max( J 4)|o = A(xq) = A max , x e ft, and is, moreover, 
not a critical point, so that VA| x=Xo ^ 0, where the gra- 
dient is taken in M. Therefore, we have that the the 



neutral manifold and the unstable manifold are tangent 
to the manifold A{x) — A max in x = xq. We have that 
the intersection between the manifolds A(x) = A and f2 is 
the empty set if A > A max . We define as S^ax the subset 
of R d included between the manifolds A(x) = A max and 
A(x) = T. Furthermore, we define as the subset 

of M. d included between the hyperplane (3 max tangent to 
the manifold A{x) = A rnax in x — xq and the hyperplane 
/3t, which is obtained by applying the translation given 
by the vector g(T— A 

max) to the hyperplane ftrnaxi 

where 
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VA\ 

X = X() 

/\VA\ 

X — Xq I * 



As T ^ A rnax , which is the 



limit of our interest, we have that Q^i = fl fl 



and fi^t 



become undistinguishable up 
to leading order. See Fig. 1 for a cartoon depicting this 
geometrical construction. 

More in general, we denote as Q v , with V > U > T, 
the intersection between and the subset of R d included 
between the hyperplane f3jj and fiy > where fix , X = U, V 
is obtained from firnax 

by applying to it the translation 
given by the vector g(X - A max ). 

It is now clear that we observe an exceedence of the 
observable A(x) above T each time the systems visits a 
point belonging to Cl^ In more intuitive terms, an 
excedeence is realized each time the system visit a point 
x € Q whose distance dist(x, (3 m ax) from the hyperplane 
Pmax is smaller than y max = (A max - T)/\WA\ x=Xo \. 

We define the exceedances for the points x <G £l T A 
as z = A(x) — T. An exceedance z corresponds geo- 
metrically to a distance y = dist(x, [3t) = z/\VA\ x — Xo \ 
from !3 T and a distance k = dist(x, /3 max ) = y max - y = 
(A max - T)/\VA\ X=X0 \ - z/\VA\ x=Xo \ from (3 max . 

Therefore, P(z > Z\z > 0) = P(y > Y\y > 0), where 
Y = Z/\VA\ x=Xo \. We have that P(y > Y\y > 0) = 
P(y > Y)/P(y > 0). In terms of the invariant mea- 
sure of the system, we have that the probability Ht{Z) 
of observing an exceedance of at least Z given that an 
exceedence occurs is given by: 



H T (Z) 



"(n£f.) 



(5) 



where we have used the ergodicity of the system. Obvi- 
ously, the value of the previous expression is 1 if Z = 0. 
The expression contained in Eq. (5) monotonically de- 
creases with Z (as ^+f x 2 C tt T A + m Z a l if Z x < Z 2 ) and 
vanishes when Z = A max — T. 



B. Derivation of the Generalised Pareto 
Distribution Parameters for the Extremes of a 
Physical Observable 

We now wish to understand how to express the nu- 
merator and denominator of Eq. 5 as a function of T, Z, 
and A max . We follow some of the ideas presented in Hol- 
land et al. [30] and use the fact that we are considering 
Axiom A systems. We define D as the local dimension 
around xq, such that the mass of the intersection of the 
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FIG. 1. A low-dimensional cartoon of the geometrical con- 
struction used for deriving the EVL for exceedances above 
the threshold T for the observable A(x) such that max(A)|n = 
A max is realized for x — xg. a) The manifolds A(x) — A max 
and A(x) — T are depicted, together with the the attracting 
invariant set Q and the two hyperplanes f3 max and /3t- fimax 
is tangent to A(x) = A max in xo and /3t is obtained from 
Pma X via translation along g(T — A max ). b) The hyperplanes 
Pm ax and Pt delimit the region E, 
Q. is S!^„ c) The manifolds A(x) 



max- Its intersection with 
A max and A(x) = T 



delimit the region E„ 



T -> A ma3: , we have that 0. max ^ma 



Its intersection with Q is fi„ 

T 



As 



ball B t (xq) and the attractor scales as r D . In the case 
of Axiom A systems, since the invariant measure is SRB 
we can write D = dn almost everywhere on the attrac- 
tor [45, 46], where dn is the Hausdorff dimension. As 
discussed in [45] , we have that in this case all definitions 
of dimension agree. Moreover, we can conjecture that 
D = dxY- where d^Y is the Kaplan- Yorke dimension 
[471: 



IKY 



AnJ-1 



(6) 



where the Aj's are the Lyapunov exponents of the sys- 
tems, ordered from the largest to the smallest, n is such 
that J2k=i is positive and J2k=i A fc is negative. Fol- 
lowing [45] , we can also write dKY — d u + d n + d s , where 
d s , d u and d n are the dimensions of the attractor f2 re- 
stricted to the stable, unstable and neutral directions, 
respectively, at the point x — Xq. We have that d u is 
equal to the number of positive Lyapunov exponents , 
d u = #({Aj > 0}, j = 1, . . . , d), d n for Axiom A systems 
is unitary, while d s is given by d s = dxY — d u — d n . Note 
that we can also express dxY as follows 



IKY 



[d s 



|A„+i| 



d u + d n + d s , (7) 



so that {d s } — Y^^ =1 Afc/|An + i|, because the last term 
gives a positive contribution smaller than 1, and d u and 
d n are both integer. 

Let's now estimate ^(^^ ) as a function of y max in 
the case of generic quadratic tangency between the the 
hyperplane A(x) — A max and the unstable manifold in 
x = xq We follow the construction proposed by Holland 
et at [30] for low dimensional maps. We derive the result 
by considering the following heuristic argument. Near xq 1 
the attractor can be seen as the cartesian product of a 
multidimensional paraboloid (of dimension d u + d n and 
of a fractal set of dimension d s immersed in R d ~ du ~ dn . 
The mass of the paraboloid scales with r du+drl , where r is 
the distance from the minimum. Instead, for each point 
of the paraboloid, the mass of corresponding fractal set 
scales with l ds , where I is the distance along the the carte- 
sian projection. In our case, I oc y m ax and r oc y/y m ax ■ 
for the former relation we assume a generic relation be- 
tween the direction of the gradient of A and the stable 
directions, while the latter relation results from the func- 
tional form of the paraboloid, see also Fig. 1). Hence, 
we obtain that the mass of il^ can be written as 
v{^A ma J oc y s max where S = d s + (d u + d„)/2. The 
argument can be made more formal by considering the 
disintegration of the SRB measure v along the stable and 
unstable directions of the flow [47] . 

As a side note, we emphasize that in the case of more 
general tangencies between the the unstable and neutral 
manifold and the manifold A(x) — A max , we have that 
5 = d s + £ J=1 d u {2j)/{2j) + d n (2j)/(2j), where d u (2j) 
(d n (2j) ) gives the number of directions along the unsta- 
ble (neutral) manifold where a tangency of order 2j with 
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the manifold A(x) = A max is found. We obviously have 
that Y,j=i du{2j) = d u . 

We continue our discussion considering the case of 
generic tangency. Since y max = (A max - T)/\VA x = Xo \, 
we have that ^(r2^ max ) oc (A max — T) s . Following the 
same argument as above, we have that v{Vi^ z ) oc 



v(n T A +z ) ex (A, 



y) , where y = Zj /|VAr =Xo | , so that we derive 
- T - Z) 5 . From Eq. 5 we obtain: 



H T {Z) 



{A r 



T - Z) s 



max 

-Tf 



A — T 

■^max 



■ (8) 



Note that the corresponding cdf is given by Ft(Z) = 
1 — Ht(Z). Comparing Eqs. 3 and Eqs. 8, one obtains 
that Ft(Z) belongs to the GPD family, and that the 
GPD parameters can be expressed as follows: 



-1/5 



(A n 



T)./5. 



(9) 
(10) 



These results complement and extend what obtained in 
[30] using the GEV approach. 



C. Comments 

Equation 9 provides a very valuable information, as 
it shows that the shape parameter £ of the GPD does 
not depend on the considered observable, but only on 
the dimensions of the stable and of the unstable mani- 
fold. Moreover, the shape parameter is always negative, 
which is hardly surprising as we are considering com- 
pacts attractor and a well-behaved observable, whose val- 
ues on the attractor have an upper bound. Note that 
for Axiom A systems, d s and d u are constant through- 
out the attractor f2, so that the information gathered 
for x — xo has a global value. Therefore, the ex- 
pression for £ is universal, in the sense that we can 
gather fundamental properties of the dynamical system 
by looking at the shape parameters of the extremes of a 
generic observables with the properties described above. 
Note also that 5 can be used to put upper and lower 
bounds to the Kaplan- Yorke dimension of the system, as 
dxv = d s + d u + d n > d s + (d u + d n )/2 = 5 = 
and d KY = d s + d u + d n < 2d s + d u + d n = 25 = -2/£, 
so that — 1/£ < dKY < — 2/£. On the other hand, these 
inequalities can be read as constraints to the shape pa- 
rameters of the extremes of a general observable for a 
system for which we know the Kaplan- Yorke dimension: 
d KY /2 < -l/£ < d KY . 

On the other hand, the expression we obtain for a pro- 
vides clear support for calling it the scale parameter. We 
derive, as anticipated, that a > and we observe that it 
is proportional to the actual range of values considered 
as extremes of the observable A(x), by incorporating the 
difference between the absolute maximum of the observ- 
able A max and the selected threshold T. Therefore, if we 
consider as observable A\(x) = aA(x), with a > and 



take as threshold for A\(x) the value aT, we have that 
£ai = £a and oa x — oca- In physical terms, a changes if 
we change the unit of measure of the observable, whereas 
£ does not. 

Let's define min( J 4)|o = A min . If we select an ob- 
servable A 2 (x) = $>(A(x)), with max($)|[ ylm . ni ^ max ] = 
$(4 raol ), $ diffcrcntiablc and d$(y)/dy positive in a 
sufficiently wide neighbourhood around y = A(x max ) 
so to ensure monotonicity of A 2 (x) near x — x max , 
we get £ A2 = £ A and a A2 = ja A , where 7 = 
d*(y)/dy\ 

y=A(x max )- 



D. From the Extremes to the Partial Dimensions 
along the Stable and Unstable Directions of the Flow 

It is worth considering the following strategy of in- 
vestigation of the local properties of the invariant mea- 
sure near x — x , where A(x ) = A max . By perform- 
ing statistical inference of the extremes of A we can de- 
duce as a result of the data fitting the best estimate of 
= 1/5. If, following [12], we select as observable, e.g. 
B(x) = — (dist(x, xq))P , P > 0, we have that the ex- 
tremes of the observable B feature as shape parameter 
£b = —P/D = — /9/dji-y[12]. We can then easily derive: 



Ka Kb 



(11) 
(12) 



where, as discussed above, we can take d n = 1. There- 
fore, using rather general classes of observables, we are 
able to deduce the partial dimensions along the stable 
and unstable manifolds, just by looking at the properties 
of extremes related to x = x n . It is to be noted that, as 
clear from the results presented in [12], similar conclu- 
sions can be drawn choosing powers of dist(x,xo) forms 
for B are possible. 



E. Expressing the shape parameter in terms of the 
GPD moments and of the invariant measure of the 
system 

We consider the physical observable A. We denote by 



Igpd{z; £ a , <ta) = -7- (Fgpd(z; £a,<7a)) 
az 



1 f 1+ Uz 



&A \ &A 



(13) 



the density corresponding to the cumulative distribution 
given in Eq. 3. We can express its first two moments as 
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follows: 

r—VAltA 

\ dz z f G PD(z;£,A,VA) 

JO 1 - C.A 



(14) 



dz z .fGPD{z;£,A,v a) 



(1-60(1-260 



= A*2- 

(15) 



It is easy to derive that 
1 



1 



Mi 



(16) 



where we indicate explicitly that we refer to the observ- 
able A. The previous equation is true both for the un- 
perturbed and the perturbed dynamics, so that: 



£a = 



Mf 2 



M 2 - Mi 



(17) 



where we have added the explicit dependence of the re- 
alized statistics on the e perturbation parameter, and we 
have introduced in the last expression the index of disper- 
sion id e A , the ratio between the variance and the squared 
first moment of the considered stochastic variable. 

We now try to connect the previous formulas to the 
properties of the invariant measure of the dynamical sys- 
tem. As we know, the GPD is the exact asymptotic 
model for the extremes of the observable A, so that we 
can express the results in terms of the conditional invari- 
ant measure as follows: 



Mi 



Jp € (dx)&(A(x)-T)(A(x)-T) _ (Atf 



J Pe (dx)e(A(x)-T) 



(A)) 



T.c 



(18) 



T , _ J p e (dx)Q(A(x) - T)(A(x) - Tf (A 2 ) T * 

A*2 — ~~ 



Jp e (dx)@(A(x)-T) 



( A oY j 



(19) 



where, in general, 

(I„) T ' e - J p e (dx)e(A(x) - T)(A(x) - T) n . (20) 

We then obtain the following expression for the shape 
parameter: 




(A 2 )^(A )^ - ((A^y 



(21) 



where the result is exact in the limit for T — > A max . It is 
useful to verify that the right hand side of Eq. 21 gives 
the same general result as given in Eq. 9. By definition 
we have: 



(Ac) 



J p e (dx)Q(A(x)-T) 
~ C T ' e (A e 



Tf 



(22) 



where the dimension 5 depends on e, while the metric 
factor £ possibly also depends weakly on e and T. We 
derive the following expression using repeatedly the dis- 
tributional relation x d/dx[Q(x)] = 0: 



so that 



<^> T ' e = -d^i> T ' e 



(23) 
(24) 

(25) 



and 



2C 



(^ + l)(<5 e + 2) 



(A max -Tf+ 2 , (26) 



where we have assumed that the dependence of ( on T 
is very weak in the regime of interest, By plugging these 
expression into Eq. 21, we indeed obtain £ e = — l/5 e , 
which agrees with Eq. 9. 

We also note that it is possible to generalize the results 
given in Eqs. 23-26. Using the fundamental theorem of 
calculus, it is possible to derive that: 



(An) 



T.< 



-f 

J z 



dzn(z-T) n - 1 (A ) z ' € 



Moreover, it is notable that what presented in this sub- 
section can be replicated step by step for the distance ob- 
servables observable B{x,Xq) = — dist(x, x^Y discussed 
above. We obtain: 




((Bi) T ' e ) 2 



2 \ (B 2 )^(A )T,* - «B 1 )r,e)2 



(27) 



where the quantities (Bj) T - e , j = 0, 1, 2 are constructed 
analogously to how described in Eq. 20. 



III. RESPONSE THEORY FOR THE 
EXTREMES OF GENERAL OBSERVABLES 

We wish to present some ideas on how to use response 
theory and the specific expressions given in Eqs 9-10 to 
derive a response theory for extremes of physical and dis- 
tance observables in Axiom A dynamical systems. Let's 
assume that we alter the Axiom A dynamical system un- 
der consideration as x — G(x) — > x = G(x) + eX(x), 
where e is a small parameter, so that the evolution oper- 
ator is transformed as /* — > /* and the invariant measure 
is transformed as p n — > Pe . Ruelle's response theory al- 
lows to express the change in the expectation value of 
a general measurable observable vE'(x) as a perturbative 
series as (*) e = (*) + YlJLi e H^^)o, with j indicating 
the order of perturbative expansion, where 



(*) e - J Pe(dx)^{x) 
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is the expectation value of over the perturbed invariant 
measure and 



Wo = J Po(dx)*(x) 

defines the unperturbed expectation value of <f\ The 
term corresponding to the perturbative order of expan- 
sion j is given by where can be expressed 
in terms of the time-integral of a suitably defined Green 
function [33]. At this stage, we limit ourselves to the 
linear response of the system. We consider the following 
useful formula: 



d n <*) ( 



de™ 



n!(* (n) > . 



e=0 



and take into account the n = 1 case. 



A. Sensitivity of shape parameter as determined 
by the changes in the moments of the distribution 

We wish to propose a linear response formula for the 
parameter £4 using Eqs. 18-21. We start by considering 
that in Eq. 21 the shape parameter is expressed for every 
T < A max as a function of actual observables of the sys- 
tem. Unfortunately, in order to apply Ruelle's response 
theory, we need the observables to be smooth, which is 
in contrast with the presence of the in the definition of 
the terms (Aj) T ' e . Nonetheless, replacing the 8's with a 
smooth approximation 0g , the Ruelle response theory 
can be rigorously applied. We now consider a sequence 
of approximating 0™ such that the measure of the sup- 
port of — 0™ is smaller than S m = (A max — T)/m. 
It is reasonable to expect that as 5 m — > 0, the effect 
of the smoothing becomes negligible, because a smaller 
and smaller portion of the extremes is affected, and the 
response of the smoothed observable approaches that of 
(Aj) T ' e . Therefore, we can retain the in the definition 
of the (Aj} T ' e and define rigorously for every T < A max : 



de 




e= — ^(A 2 )^ {Ao) T,e_ {{Al) T,e ) 2 ) ^ 

(28) 

By expanding the derivative in Eq. 28, the previous ex- 
pression can be decomposed in various contributions en- 
tailing the linear response of the system to the e pertur- 
bation for the observables (A ) T , (Ai) T , (A 2 ) T and their 
values in the unperturbed case for e = 0. 

Therefore, Eq. 28 provides a recipe for computing the 
sensitivity of £ A ' 
terest, where A. 
to collect experimental data or process the output of nu- 
merical simulations we need to select a threshold which 
is high enough for discriminating true extremes and low 
enough for allowing a sufficient number of samples to be 
collected for robust data processing. Note that all statis- 
tical procedures used in estimating GPD parameters from 



at e = for any case of practical in- 
max — T is indeed finite, because in order 



data are actually based on finding a reasonable value for 
T such that both conditions described above apply by 
testing that parameters' estimates do not vary apprecia- 
bly when changing T [4]. So, if when investigating the 
extremes of A for the unperturbed and e— perturbed dy- 
namics we find a common value of T such that the GPD 
statistical inference of extremes is satisfactory, Eq. 28 
provides the correct formula for the sensitivity of £^' c at 
e = 0. 

We wish to underline that apparently formal problems 
emerge when taking the limit if Eq.28 for higher and 
higher values of T. It is indeed not clear at this stage 
whether 



lim 

}A ma ^ 



de 



e=0 



= lim 

T^A mt 



f 1 '' 



C T,0 



, e->0 



(29) 



exists, because we cannot apply the smoothing argument 
presented above in the limit of vanishing A max — T . More- 
over, it is not clear whether such limit is equal to 



lim lim 

e->0 T — yA m( 



C T,0 
$A 



(30) 



which seems at least as well suited for describing the 
change of the shape observable given in Eq. 21 due to an 
e— perturbation in the dynamics. Obviously, if the two 
limits given in Eqs. 29 and 30 exist and are equal, then 
a rigorous response theory for £4 can be established. 

The same derivation and discussion can be repeated 
for the B observables introduced above. In this case, we 
can define:: 



de 




«£i> T < e ) 2 



£=0 _ ^ 1J2r ^B )T,e- {{Bl )T,e ) 2 ) ^ 

(31) 

where the correct limit for T is T — > 0. 

Let's try to give a more intuitive interpretation to the 
results given above. Let's consider Eq. 17 and assume 
that, indeed, £ is differentiable with respect to e. We 
have: 



de 



Id 
2ch 



1 

id*-. 



1 



. -0 - >il 1 1 ) -0 2id e A 2 de 

(32) 

which implies that the sensitivity of the shape parameter 
is half of the opposite of the sensitivity of the inverse of 
the index of dispersion id a- Therefore, a positive sensi- 
tivity of the index of dispersion (larger relative variability 
of the extremes of the observable A with positive values 
of e) implies a larger value (closer to 0) of £4, and so the 
possibility that larger and larger extremes are realized. 
Same interpretation applies for the B observables. 



{id A } 



B. Sensitivity of the shape parameter as 
determined by the modification of the geometry 

In the previous subsection we have shown that the Ru- 
elle response theory supports the idea that the shape pa- 
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rameters descriptive of the extremes of both the physi- 
cal observables A and the distance observables B change 
with a certain degree of regularity when considering 
e— perturbations to the dynamics. 

In this subsection, we wish to look at the sensitivity 
of extremes with respect to perturbation from another 
angle, i.e. through the relationship between the shape 
parameters £4 and £ B and the partial dimension of the at- 
tractor along the stable, neutral and unstable manifolds 
of the underlying dynamical system, see Eqs. 11-12. As 
long as the e-perturbation is small, the modified dynam- 
ical system belongs to the Axiom A family, so that the 
results presented above apply. Therefore, we can write 
in more general terms: 

^ = -1/^ = -1/(^ + ^/2 + ^/2) (33) 
& = -0/d*y = -0/(<5 + + (34) 

In the following, we introduce somewhat carelessly 
derivatives with respect to e of quantities that are not, a 
priori, differentiable. The main point we want to make 
is that if £a and £b are differentiable with respect to e, 
then various quantities describing the structure of the at- 
tractor are also differentiable. Therefore, the existence of 
the limits given in Eqs. 29 and 30 (and their equivalent 
for the B observables) would have far-reaching conse- 
quences. We will discuss the obtained results at the end 
of the calculations. Another caveat we need to mention 
is that Eqs. 33-34 are in general true almost anywhere, 
so that we may have to interpret the derivatives in this 
section in some suitable weak form. 

It seems relevant to add the additional hypothesis of 
strong transversality for the unperturbed flow, which is 
equivalent to invoking structural stability [47]. We take 
such pragmatic point of view and proceed assuming that 
derivatives with respect to e are well defined. Linearizing 
the dependence of £a on e around e = in Eq. 33, we 
obtain: 



de 



d(dt) de 



£=0 



d(d«) de 
f dq d(Q 
\d(d*) de 



£=0 



£=0 

(35) 



We have that d(c^)/de| e= o = d(c^)/de| e= o = 0, as, 
thanks to structural stability, small perturbations do not 
alter the qualitative properties of the dynamics, and can- 
not change in a step-wise way the number of expanding 
or neutral directions. We now separate the quantity d s 
into its integer component and the rest: 



d^ 
de 



= ( d& d«) 
U(di) 



de 



d& d(K» 
d(d«) de 



£=0 

Only the last term is different from zero, because, thanks 
to structural stability, d([d*])/de| e= o = 0. Using Eq. 7, 



we obtain: 



da 

de 



de 



K+il 



£=0 



d(dj) de 



(37) 



£=0 



where n is defined as in Eq. 6. Since n is not altered for 
infinitesimal e perturbations to the dynamical system, we 
have that, using Eq. 6: 



a \K+i\ 



de 



d(d* KY ) 



de 



(38) 



Expanding the previous expressions, the final formula 
reads as follows: 



da 

de 



£=0 



1 



d{d* KY ) 



(d«+d«/2 + d*/2) 



de 



(39) 



£=0 



This implies that the shape parameter £ increases, thus 
attaining a value closer to zero (£4 is always negative) 
when the perturbation increases the Kaplan- Yorke di- 
mension of the attractor, so, in qualitative sense, if it 
favors "forcing" over "dissipation". This matches quite 
well, at least qualitatively, with the discussion following 
Eq. 32. 

We have that, when considering a distance observable 
of the form B{x) = — dist(x, xqY , following the same 
steps described above one gets the following result: 



de 



= ( P d(rf^) \ 
e= o \d e KY 2 de J 



(40) 



e=0 



such result can be easily generalized by considering the 
class of observables described in [12]. 

Combining Eq. 28 with Eq. 39 and Eq. 31 with Eq. 
40, we can derive two expressions for the derivative of 
the sensitivity of the Kaplan Yorke dimension at e = 0: 



d(d* KY ) 



de 



£=0 



(dl + di/2 + dl/2f 



£=0 



«4 



de (A 2 )^(A )^ - ({A 1 ) T ' e ) 2 



£=0 
(41) 



KY 



£=0 



d ((B^) 2 
de (b 2 )t, £ (b )t, £ _ ((B 1 )T,e)2 



£=0 
(42) 



where we take the limit for T — y A max in Eq. 41 and 
T -> in Eq. 42. 

The previous results imply that if one of £b or 
the Kaplan- Yorke dimension of the underlying Axiom A 
system change smoothly with e— perturbations to the dy- 
namics, so do the other two quantities. This may sug- 
gest ways to study the regularity of the Kaplan- Yorke 
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dimension by resorting to the analysis of the regularity 
of a much simpler expressions involving moments of given 
observables. 

This result provides useful insight also not considering 
the problematic limits discussed above. Taking a more 
qualitative point of view, It suggests that when consider- 
ing small perturbation in the dynamics of chaotic systems 
behaving like Axiom A systems there is a link between 
the presence (or lack) of regularity of the parameters de- 
scribing the extremes of a wide class of observables and 
of the regularity of the Kaplan- Yorke dimension. This t 
matches with the fact that numerical practice with mod- 
erate to high dimensional strongly chaotic systems shows 
that, actually, the parameters describing extremes of en- 
ergy like observables and the Lyapunov exponents both 
have - within numerical precision - a smooth dependence 
on the parameters of the system. See [3, 48] for an exten- 
sive discussion in a simplified yet relevant fluid dynamical 
model. A detailed investigation of the apparent regu- 
larity for all practical purposes of the Lyapunov expo- 
nents with respect to small perturbations in the dynam- 
ics of intermediate complexity to high-dimensional mod- 
els has been presented in [49]. We also wish to remark 
that if these regularity hypotheses were not satisfied, the 
very widespread (and practically successful) procedure 
of parametric tuning of high-dimensional models of nat- 
ural, engineered or social phenomena would be absolutely 
hopeless, and delicate numerical procedures such as those 
involved in data assimilation of high-dimensional dynam- 
ical systems would lack any sort of robustness, contrary 
to the accumulated experience. 



IV. CONCLUSIONS 

This paper has addressed the problem of studying the 
EVT for general observables of mixing Axiom A dynam- 
ical system. We have set in a common framework the 
investigation of general properties of distance observables 
B, for which we had derived some basic results in [12], 
and of physical observables, first discussed in [30]. By 
physical observables we mean rather general classes of 
smooth observables A (e.g. quadratic, energy-like quan- 
tities) which achieve an absolute maximum A max at a 
specific, non-critical (the gradient does not vanish) point 
of the attractor x = xq of the system. 

We have built up from the recent results of Holland 
et al. [30], who have studied accurately this problem on 
discrete maps using the GEV framework and have de- 
rived, as fundamental result, that the extremes of A in- 
deed obey a GEV extreme value law and that £4, the 
shape parameter of the distribution, is equal to —1/5, 
with 5 = d s + d u /2, where d s is the partial dimension 
along the stable manifold and d u is the partial dimension 
along the unstable manifold. Using the GDP approach, 
thus considering exceedances above a given threshold T, 
and considering the physically relevant case of mixing 
Axiom A systems, we derive through direct integration 



that the shape parameter £a can also be expressed as 
£a = — 1/5- We have framed our results for continuous 
flows, so 5 = d s +d u /2 + d n /2, where d n is the dimension 
along the neutral direction and is unitary. We need to re- 
mark that such results has been constructed using some 
intuitive geometrical construction and assuming generic 
relation between the direction of the gradient of A for 
x — x and the direction of the stable directions. While 
it is possible to construct fractal sets where such generic 
conditions are not realized, we assume that typical Ax- 
iom A (or more loosely, sufficiently chaotic systems) obey 
instead such conditions. This can be tested numerically. 

On the other hand, we have been able to derive the 
explicit expression for the scale parameter a = (A max — 
T)/5. It is clear that the £ parameter is always negative 
(so that the distribution of extremes is upper limited), 
reflecting the fact that the observable is smooth and the 
attractor is a compact set. Moreover, measuring £ allows 
is to provide an upper and lower bound for dxy and vice 
versa, because —1/£a < d KY < —2/^, or, conversely 
d KY /2 < -l/U < d KY - 

In particular, we have that £a is small and negative 
if and only if very the Kaplan- Yorke dimension of the 
attractor is large. If we consider a system with a high 
dimensional chaotic system (e.g. in the case of an ex- 
tensive chaotic system with many degrees of freedom), 
we derive that £4 ~ 0. This may well explain why in 
a multitude of applications in natural sciences such as 
hydrology, meteorology, oceanography the special £ = 
member of the GPD - the exponential model - given in 
Eq. 3 usually gives a good first guess of the statistics of 
observed extremes [9]. 

Alternatively, this result suggests that if we perform 
a statistical analysis using the POT method (using 
an empirical threshold T) of the extremes for a high- 
dimensional chaotic system and obtain as a result of the 
statistical inference of the collected data for the GPD 
model a shape parameter £4 <C or £4 > 0, we should 
conclude that our sample is not yet suited for an EVT 
statistical fit. This may depend on the fact that we have 
selected an insufficiently stringent value for T. Obviously 
choosing higher values for T implies that we need to have 
longer time series of the observable under investigations. 

Interestingly, by combining the expression for £a ob- 
tained in this paper for a physical observable A and the 
expression for £g of the GPD describing the extremes 
of observable of the form B(x) — — dist(x, xo) , where 
dist is the distance function [12], we can express the par- 
tial dimension of the attractor on the stable and unsta- 
ble dimension as simple functions of £4 and ^b- The 
straightforward results is that 2/£a — 2/3/£b = d u + d n 
and /?/£b — 2/£a = d s - This provides further support 
to the idea that extremes can be used as excellent diag- 
nostic indicators for the detailed dynamical properties of 
a system. The message seems to be that one can con- 
struct observables whose large fluctuations give precise 
information on the dynamics. 

In the second part of the paper we have tackled the 
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problem of studying how the properties of extremes of 
the observable A change when a e-perturbation is added 
to the system. As theoretical framework wc have taken 
the point of view of Ruelle [33, 34], who has shown that 
the SRB measure of Axiom A systems is differentiable 
with respect to e-perturbations to the dynamics and has 
provided explicit formulas for studying how the expec- 
tation value of generic observables of Axiom A systems 
change when the system is subjected to perturbations. 

We have used the fact that the GPD is an exact asymp- 
totic model for extreme events of the observable A to find 
a simple functional relation between £4 and £g and the 
first two moments of the probability of above-threshold 
exceedence expressed in terms of the invariant measure 
of the system perturbed with an when the system for a 
general perturbation The result we obtain confirms the 
output of the direct calculation on the scaling proper- 
ties of the measure near x - By expressing £4 and ^in 
terms of a simple combination of moments of the em- 
pirical distributions of above-threshold exceedances, we 
find expressions which are amenable to direct treatment 
with Ruelle's response theory, at least when we do not 
consider the limit T — > A max but stick to the practical 
situation where we need to consider a finite range for the 
extremes. The differentiability properties of xia and xis 
are hard to ascertain in the limit, 

We have then taken into consideration our results on 
the relationship between £4 and £b and the partial di- 
mensions of the attractor. Interestingly, it seems that 
there is an intimate connection between the differentia- 
bility with respect to e of £a, £b and dxy, so that ei- 
ther all of them or none of them is differentiable with 
respect to e. Under the hypothesis of differentiability, we 
have been able to derive that the sensitivity of £4 and 
£b with respect to e is proportional to the sensitivity of 
the Kaplan- York dimension (Iky to e, in such a way that 
if the perturbation tends to increase the dimensionality 
of the attractor (thus, in physical terms, favoring forcing 
over dissipation), the value of £ becomes closer to zero, so 
that the occurrence of very large extreme events becomes 
more likely. The system, in this case, has more freedom 
to perform large fluctuations. 

Taking a more pragmatic point of view, these results 
at least provide a rationale for the well-known fact that 
in moderate to high-dimensional strogle chaotic systems 
the Kaplan- Yorke dimension (and, actually, all the Lya- 
punov exponents) change smoothly with the intensity of 
the perturbation vector field, as discussed in [48, 49], and 
similar behavior is found for the parameters describing 
the extremes of energy-like quantities [3] . 

Since it has been shown in [36] that the linear response 
vanishes for any observable in the case of stochastic forc- 
ing of rather general nature, whereas the second order 
response gives the leading order of perturbation, we ex- 



pect that adding a moderate noise to a chaotic dynamical 
system will not alter significantly the shape parameter £ 
describing the EVT of both physical and distance observ- 
ables. 

In this work we have considered the case where the 
observable A has a unique maximum restricted to f2 in 
x = x e £1. If ft and A share some symmetries, x 
is not unique, and instead there is a set of points x 's 
belonging to f2, finite or infinite, depending of the kind 
of symmetries involved, where A reaches its maximum 
value restricted to O. Let's consider the relevant case 
that A and fi share a discrete symmetry, so that 
the set of the maximal point xo's, has finite cardinal- 
ity. The results discussed here for the extremes of A will 
nonetheless apply, because we can perform an equivalent 
geometrical construction as in Fig. 1 for each element 
of xo- When we consider an e-perturbation to the dy- 
namics which respects the discrete symmetry, it is clear 
that all the results of the response presented here apply. 
Finally, one can deduce that if the considered perturba- 
tion, instead, breaks the discrete symmetry, the results 
presented here will still be valid as the break of the de- 
generacy will make sure that only one of the xo's (or 
a subset of xo ; if the corresponding the perturbed vec- 
tor flow obeys to a a subgroup of the original symmetry 
group) is still accounts for the extreme events of A. 

This works may constitute a new theoretical viewpoint 
for studying extremes in a rather general setting and un- 
derstanding how they change when the dynamical system 
is slightly perturbed. The results on how to express the 
sensitivity of the Kaplan- Yorke dimension with respect 
to an e— perturbation to the dynamics seems also useful. 
Our findings have relevance in a multitude of sectors of 
physical, engineering and social sciences, and, just to pro- 
vide a basic example in the investigation of the impact of 
climate change on climate extremes. Obviously, as in the 
case of all results pertaining to EVT, it is important to 
test numerically the practical verification of the findings 
presented here. 
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